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1. Introduction 

The muon anomalous magnetic moment g — 2 [a M = (g — 2)/2] has been measured 
by the E821 experiment (Muon g-2 Collaboration) at Brookhaven National Labo- 
ratory (BNL) with an impressive accuracy of 0.72 ppm^ yielding the present world 
average^ 

a c * p = 11 659 208.0(6.3) x lO" 10 (1) 

with an accuracy of 0.54 ppm. New experiments^ 1 are under design with a goal of 
measuring a M with an accuracy of at least 0.25 ppm. 

On the theory side, a large amount of work has been devoted to reduce the 
uncertainty of the Standard Model prediction. A recent updated discussion with an 
extensive list of references for both theoretical predictions and experimental results 

*On leave of absence from CAFPE and Departamento de Fi'sica Teorica y del Cosmos, Universidad 
de Granada, Campus de Fuente Nueva, E-18002 Granada, Spain. 
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is Ref. 3 and Ref. [H and a more introductory exposition can be found in the lectures 
by Knecht in Ref. O 

In this paper, we review the present status of the hadronic light-by-light con- 
tribution (hLBL). A somewhat shorter version is the published talk in Ref . El The 
uncertainty in the hLBL is expected to eventually become the largest theoretical 
error. This contribution is shown schematically in Fig. [I] It consists of three photon 
legs coming from the muon line connected to the external electromagnetic field by 
hadronic processes. Its contribution can be written as 




P P4 P5 P 



Fig. 1. The hadronic light- by-light contribution to the muon g — 2. 
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q 2 p\p\{p\ — to 2 ) (p\ — to 2 ) 
x W v ^{p l ,p^p i )u{p , ) la {^ + m) lv {^ +m) lp u[p) (2) 

where q — p\ +P2 +P3- To obtain the amplitude M. in ([2]), the hadronic contribution 
to the full correlator n pra/3 {p\ , p2 , P3 — > 0) needs to be known for all possible four- 
momenta pi and P2 ■ The correlator is defined via 

n pva P(pi,p2,p 3 ) = i 3 J d 4 x J d 4 y J d A ze i{pl - x+P2 - y+ps - z) x 

<0|T [^( )^(x)^(t / )^ 3 (^)] |0> (3) 

with V*(x) = qQj^q (x) and Q = diag(2, — 1, — 1)/3 the quark charges. The 
external magnetic field couples to the photon leg with momentum P3 — > 0. In the 
remainder whenever we refer to we specifically mean only the hadronic light-by- 
light contribution to it. 

Clearly, the correlator ([3]) is a complicated object. It contains many independent 
Lorentz structures, each of comes with a function of the variables p\, p\ and q 2 . As 
a consequence, many different energy scales can be involved in the calculation of the 
hadronic light-by- light contribution to muon g — 2. This makes it difficult to obtain 
the full needed behavior of the correlator (El) from known constraints. Therefore 
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no full first principles calculation exists at present. The needed results cannot be 
directly related to measurable quantities either. Lattice QCD calculations are at 
the exploratory stage only, see e.g. Ref. [7] 

In fact, there has long been a confusion about hadronic exchange^ versus quark 
loop estimates. This confusion was resolved by organizing the different contributions 
according to the lowest power in 1/N C and the lowest order in the chiral perturba- 
tion theory (CHPT) expansion counting where they start contributing^. One can 
distinguish four types of contributions: 

• Goldstone boson exchange contributions are order N c and start contributing 
at order p e in CHPT. 

• (Constituent) quark-loop and non-Goldstone boson exchange contributions 
are order N c and start contributing at order p & in CHPT. 

• Goldstone boson loop contributions are order one in 1/N C and start con- 
tributing at order p 4 in CHPT. 

• Non-Goldstone boson loop contributions are order one in 1/N C and start to 
contribute at order p 8 in CHPT. 

The two existing full calculations^^, are based on this classification. The Gold- 
stone boson exchange contribution (GBE) was shown to be numerically dominant 
in Refs. [8] and [TO] after strong cancellations between the other contributions. But 
the other contributions, though each smaller than the GBE, were not separately 
negligible. Using effective field theory techniques, Ref. ITT1 showed that the leading 
double logarithm comes from the GBE and was positi ve. Refs. [XT] and [12] found a 
global sign mistake in the GB E of the earli er wor k^-^ which was confirmed by the 
authors of those wor and by otherPHa. In the remainder we will always 

correct for this sign mistake without explicitly mentioning it. 

Recently, Melnikov and Vainshtein pointed out new short-distance constraints on 
the correlator ([3|) in Ref. 17, studied and extended in Ref. HH The authors of Ref. [17| 
constructed a model which satisfies their main new short-distance constraints in 
order to study its effects and found a number significantly different from the earlier 
work. They approximated the full hLBL by the GBE and axial-vector exchange 
contributions. 

One of the purposes of this review is to critically compare the different contribu- 
tions in the different calculations and extend somewhat on our earlier comments^. 
For the dominant GBE we also present some new results on the momentum regions 
which are relevant. In earlier work several studies were done to check which mo- 
mentum re gion s were important. These used different methods, varying the vector 
meson masi^, studying the cut-off de pende nce^ and expansions around various 
momentum regions in the loop integral d^^ . 

In Sect. [5] we discuss the calculations done before 2002 and compare their results. 



a Wc stick here to the formulation "exchange" as used by us^l It is often referred to as "pole" 
contributions. We consider this misleading because the exchanged particle is used off-shell. 
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Fig. 2. A generic meson exchange contribution to the hadronic light-by-light part of the muon 



Sect. [3] discusses the short distance constraints proposed by Melnikov and Vainshtein 
and the numerical results presented in their paper. In Sect. [4] we show in detail in 
which momentum regions the contributions from 7r° exchange originate and Sect. [5] 
compares and comments on the various contributions in the different calculations. 
Finally, we present our conclusions as to the present best value and error for the 
hLBL in Sect. © 

2. Results obtained up to 2002 

In this section we discuss the calculations performed in the period 1995-2001. These 
were organized according to the large iV c and CHPT countings ^ discussed above. 
The CHPT counting is used as a classification tool, none of these calculations were 
actually performed at a fixed order in CHPT. We want to emphasize once more 
that the calculations in Refs. 151 ITDl H"3l and [T4l showed that only after several large 
cancellations in the rest of the contributions, the numerically dominant one is the 
Goldstone boson exchange. In this section we concentrate on the work in Refs. [8J 
and [13l with some comments and results from Refs. [101 [T2] and HH 

2.1. Pseudo- Scalar Exchange 

The pseudo-scalar exchange was saturated by the Goldstone boson exchange in 
Refs.0[IOl[l3]and[H This contribution is shown in Fig. [2] with M = n°,r],r)'. 
Refs 151 and fT3l used a variety of 7r°7*7* form factors 



fulfilling several possible QCD constraints. A more extensive analysis of this form 
factor was done in Rcf. 20 finding very similar numerical results. In particular, the 
three-point form factors used in Refs. 151 and [T31 had the correct QCD 



9-2. 



(4) 
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Table 1. Results for the ir° and it , r\ and n' exchange contributions. 

10 10 X a u 



Bijnens, Pallante and Prades 
Hayakawa and KinoshitaE3H] 
Knecht and NyffelerDSl (h 2 = 0) 
Knecht and NyffelerCS (h 2 = -10 GeV 2 ) 
Melnikov and Vainshtein H3 



7T° only 


7T°, rj and n' 


5.6 


8.5 ± 1.3 


5.7 


8.3 ± 0.6 


5.8 


8.3 ± 1.2 


6.3 




7.65 


11.4±1.0 



Table 2. Results for the axial-vector exchange contri- 
butions. 



Axial- Vector Exchange Contributions 


10 10 x a M 


Bijnens, Pallante and Prades ^ ^ 
Hayakawa and Kinoshita HSIH1 
Melnikov and Vainshtein li^J 


0.25 ± 0.10 
0.17 ± 0.10 
2.2±0.5 



short-distance behavioo 

T(Q 2 ,Q 2 )^ A/Q 2 , T{Q 2 ,0)^B/Q 2 , (5) 

when Q 2 is Euclidean. These form factors were in agreement with available data 
including the slope at the origin as well as treating the 7r°, rj and rf mixing. All form 
factors converged for a cutoff scale A ~ (2 — 4) GeV and produced small numerical 
differences when plugged into the hadronic light-by-light contribution. 

Somewhat different J-(p\,p\) form factors where used in Refs. QUI H"2l and fT4l 
but the results agree well. For comparison, one can find the results of Refs. [HI EH 
H2l IT3l and [14] in Tab. Q] of the ir° exchange and after adding rj and 77' exchange 
contributions to the dominant 7r° one. 



2.2. Axial- Vector Exchange 

This contribution is depicted in Fig. [2] with M = A = a?,/i and possibly other 
axial-vector resonances. For this contribution one needs the A77* and ^7*7* form 
factors. Little is known about these but there exist anomalous Ward identities which 
relate them to the -P77* and P^f*^* form factors. 

This contribution was not studied by Knecht and NyfFeleffl Refs.[8l [T0l[T3l ands 
[14] used nonet symmetry, which is exact in the large N c limit, for the masses of the 
axial- vector resonances. Their results are shown in Tab. [2] for comparison. 



b The observance of QCD short-distance constraints was implemented for this one and several 
other contributions in Refs. l8l and 1 131 contrary to the often heard wrong claim that Ref. 1171 is the 
first calculation to take such constraints into account, e.g. see Ref. 1211 
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Fig. 3. Quark-loop contribution, as modeled in ENJL. 

2.3. Scalar Exchange 

This contribution is shown in Fig. [2] with M = S = ao , /o and possible other 
scalar resonances. For this contribution one needs the S77* and £7*7* form factors. 
Within the extended Nambu-Jona-Lasinio (ENJL) model used in Refs. [5] and [TBI 
chiral Ward identities impose relations between the constituent quark loop and 
scalar exchanges. The needed scalar form factors are also constrained at low energies 
by CHPT. Refs. 181 and [TBI used none t symmet ry for the masses. This contribution 
was not included by the other groupi^SlHEI. 

The leading logarithms of the scalar exchange are the same as those of the pion 
exchange but with opposite sigrffSl. Refs. 151 and [TBI find that sign for the full scalar 
exchange contribution, obtaining 

a M (Scalar) = -(0.7 ± 0.2) x 10" 10 . (6) 

2.4. Other contributions at leading order in 1/N C . 

This includes any contributions that are not modeled by exchanged particles. At 
short-distance, the main one is the quark-loop. At long distances they are often 
modeled as a constituent quark-loop with form factors in the couplings to photons. 
This corresponds to the contribution shown in Fig. [31 Refs. [5] and [TBJ split up the 
quark momentum integration into two pieces by introducing an Euclidean matching 
scale A. At energies below A, the ENJL model was used to compute the quark-loop 
contribution while above A a bare (partonic) heavy quark loop of mass A was used. 
The latter part scales as 1/A 2 and mimics the high energy behavior of QCD for 
a massless quark with an IR cut-off around A -see footnote b . Adding these two 
contributions yields a stable result as can be seen in Tab. [3] 

2.5. NLO in 1/N C : Goldstone Boson Loops 

At next-to-leading order (NLO) in 1/N C , the leading contribution in the chiral 
counting to the correlator in @, corresponds to charged pion and Kaon loops 
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Table 3. Sum of the short- and 
long-distance quark loop contributionJS 
as a function of the matching scale A. 



A [GeV] 


0.7 


1.0 


2.0 


4.0 


10 10 x a M 


2.2 


2.0 


1.9 


2.0 



Table 4. Results for the charged and Kaon loop 

contributions to the hadronic light-by-light contribution to muon g — 2. 

Charged Pion and Kaon Loop Contribut ions 10 10 X a M 

Bijnens, Pallante and Pradcs (Full VMD) -1.9 ± 0.5 

Hayakawa and Kinoshita (HGS) l lu l 14 l -0.45 ± 0.85 
Melnikov and Vainshtein (Full NLO in 1/AT C guess) El Oil 

which can be depicted analogously to the quark-loop in Fig. [3] but with charged 
pions and Kaons running inside the loop instead. In general one expects loops of 
heavier particles to be suppressed and has only been evaluated for the pion loop 
and the much smaller Kaon loop. 

The needed form factorsEI r y*P + P~ and j*j*P + P~ vertices were studied exten- 
sively in Ref. [H In particular which form factors were fully compatible with chiral 
Ward identities were studied. The full vector meson dominance model (VMD) is one 
model fulfilling the known constraints. The conclusion unfortunately is that there 
is a large ambiguity in the momentum dependence starting at order p 6 in C HPT. 
Both the full VMD modePSl and the hidden gauge symmetry (HGS) modeli^JJ 
satisfy the known constraints. Unfortunately, this ambiguity cannot easily be re- 
solved since there is no data for 7*7* — » ir + ir~ . Adding the charged pion and Kaon 
loops, the results obtained in Refs. 151 and IT01 are listed in Tab.|U 

In view of this model dependence, the authors of Refs . 151 and [131 considered that 
the difference between the results from Ref. 8 and Ref. [TU] for this contribution 
needs to be added linearly to the final uncertainty of the hadronic light-by-light 
contribution to a M . 

3. New Short-Distance Constraints 

Melnikov and Vainshtein pointed outP^ a new short-distance constraint on the cor- 
relator (J3j) . This constraint is for 

(T[V^ Pl )V a (p 2 )VP(-q = -vi -pa)]|7(P3 -> 0)) (7) 



c Note that neither the ENJL model nor any fixed order in CHPT was used in any of the estimates 
of this contribution. 
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Table 5. Results quoted in Ref. 1171 for the pseudo- vector exchange 
depending of the /i(1285) and /i(1420) resonances mass mixing. 



Mass Mixing 


10 10 X a M 


No OPE and Nonet Symmetry with M=1.3 GeV 


0.3 


New OPE and Nonet Symmetry with M= 1.3 GeV 


0.7 


New OPE and Nonet Symmetry with M= M p 


2.8 


New OPE and Ideal Mixing with Experimental Masses 


2.2 ± 0.5 



and follows from the OPE for two vector currents when Pf ~ Pf 3> Q 2 with 
Pi = -Pi Pi = -Vl and Q 2 = ~q 2 : 

T[V^ Pl )V a (p 2 )} ~ E va ^ (p^/p 2 ) [qQ 2 lf3 l 5 q}( Pl + p 2 ) (8) 

with p = (pi — p 2 )/2 ~ pi ~ — p2 and Q is the light quark electrical charge matrix 
([3|) . This constraint was afterward generalized in Ref. HH Note that the new part is 
the use of {H} for the full correlator ((3]). Short-distance was already used to obtain 
the first constraint in ([5]). 

The authors of Ref. [T7| saturated the full correlator by exchanges. The new OPE 
constraint is satisfied by introducing a pseudo-scalar exchange with the vertex on 
the q,P3 side of Fig. [2] point-like rather than including a form factor. This change 
strongly breaks the symmetry between the two ends of the exchanged particle. There 
are also OPE constraints for P 2 « P 2 « Q 2 and Pf « Q 2 P 2 , essentially derived 
from the quark-loop behavior in this regimi^. Both latter OPE constraints on the 
correlator ([3]) are not satisfied by the model used in Ref. [17] but they argued that 
this made only a small numerical difference of order 0.05 x 10~ 10 . 

Ref. 17 added to the pseudo-scalar exchange an axial- vector exchange contribu- 
tion. They found this contribution to be extremely sensitive to the mixing of the 
resonances /i(1285) and /i(1420) as can be seen in Tab. [S] taken from the results 
there. The difference between the lines labeled "No OPE" and "New OPE" is the 
effect of making the q,P3 vertex point-like. The authors of Ref. Q~7]took the ideal 
mixing result for their final result for a^. 

4. Momentum Regions for tt° Exchange 

We were somewhat puzzled by the effect when saturating the new short distance 
constraint by GBE in Ref. [TT] and have therefore done a few studies to see whether 
the changes there come from large momentum regimes or are located elsewhere. 
This was because our total estimate of the quark-loop was similar to the numerical 
change in the GBE of Ref.[17l In order to do this study, we have adapted the method 
used in Refs. [8l and l20l to various form factors used in earlier works. We rotate the 
integrals in ([2]) into Euclidean space. The eight dimensional integral can be easily 
reduced to a five dimensional integral. Here one can choose as variables^ Pi , P 2 and 

d In Refs. [8l and 1201 a different set was used not quite as suitable for the present study. 
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Table 6. 7r°-exchange results for 10 10 X a M with a cut-off on the 
three photon momenta for the four cases described in the text. The 
last column is the difference between MV and KN form factors. 
The numerical error is at or below the last digit quoted. 



Cut-off A (GeV) 


WZW 


VMD 


KN 


MV 


MV-KN 


0.5 


4.74 


3.37 


3.39 


3.68 


0.29 


0.7 


7.51 


4.41 


4.47 


5.01 


0.54 


1.0 


11.3 


5.14 


5.29 


6.15 


0.86 


2.0 


21.9 


5.60 


5.99 


7.34 


1.35 


4.0 


33.8 


5.65 


6.20 


7.79 


1.59 


8.0 


49.6 


5.65 


6.24 


7.92 


1.69 


16.0 


68. 


5.64 


6.23 


7.96 


1.73 



Q and two angles 9\ and 62- These are the angles between the Euclidean p±, P2 and 
the muon momentum while Pi, P2 and Q are the size of the Euclidean momenta 
with P 2 = —p\, P 2 2 = —p 2 an( l Q 2 — —q 2 - Ref- 12 performed the integrals over three 
of these quantities analytically but not for the asymmetric case used by Ref. [17] 
We have th erefor e use numerical integration. The main integration routine used 
by us earlieJ^Q] was VEGAS. For the present study we have also performed the 
integration using an adaptive Gaussian multidimensional integration routine and 
have checked for several quantities that both agree and reproduce earlier known 
results. 

We will show the contributions to the muon anomalous magnetic moment from 
7r° exchange for several different form factors. These correspond to the point-like 
7r°7*7* form factor (WZW), the full vector meson dominance model (VMD), the 
LMD+V form factor^ with h 2 = -10 GeV 2 (KN) and the latter form factor but 
with the point-like version on the soft-photon end^ (MV). We will refer to these 
form factors as WZW, VMD, KN, and MV in the remainder of this section. We 
have used the values hi = 0, /15 = 6.93 GeV 2 and the value of /17 as given by Ref. 
IT21 We picked the value of /12 that was arguecP^ to better produce subleading OPE 
constraints. It raises the central value somewhat compared to /12 = as shown in 
Tab. Q] 

As inputs we used M v = M Vl = 0.770 GeV and My 2 = 1.465 GeV, F„ = 
92.4 MeV and the measured 7r° and muon masses. This is the origin of the minor 
differences with Ref 1 121 

As a first indication where the contributions to a M come from, we have listed in 
Tab.[6]the value of a M for the four cases with the constraint Q, Pi, P 2 < A. We have 
shown the logarithmically square divergent point-like case here to show the size of 
the suppression introduced by the form factors. Note that we cannot reproduce the 
7.65 of Ref. i 171 but we do reproduce the results of Rcfs. 8, 12landl2"01 The new short- 
distance constraint © came from the region Q ^ Pi w P 2 . We have thus checked 
how much of the difference and total comes from the region with Q < min(Pi, Pa) 
and from the region with Q larger than at least one of (Pi,P2), the numbers quoted 
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a£ P (MV) 



7e 


-09 r 


6e 


■09 - 


5e 


■09 - 


4e 


-09 - 


3e 


■09 - 


2e 


-09 - 


1e 


-09 - 




ft 




Fig. 4. The quantity a pp of Eq. 1 10 1) as a function of Pi and P2 for the MV choice. 

are for A = 16 GeV. The numbers are 10 10 x a M . 

Q < min(P 1 ,P2): 4.01 (KN) 4.74 (MV) 0.73 (MV-KN) 

Q > min(Pi,P 2 ): 2.24 (KN) 3.23 (MV) 0.99 (MV-KN) l '" 

As one sees, in fact, most of the difference comes from the region where the OPE 
condition is strongly violated. 

The results in Tab. [5] give only a partial indication of which momentum regions 
are important. In the remaining figures we therefore show the contribution to 
in several ways. We always denote p\,P2 as the momenta on the ir° side with both 
photons connected to the muon line and q the momentum on the soft-photon side. 
We can thus rewrite the contribution to a M of @ in various ways: 

a M = J rfPirfPs a p /{P 1 ,P 2 ) 
= J dhdh afi L (l u l 2 ) 

= fdhdl 2 dl q a^(l u l 2 ,l q ), 

with h = log(Pi/GeV), h = log(P 2 /GeV), and l q = log(Q/GeV) . (10) 

In Fig. [3] we have plotted 10 10 x a pp (Pi, P2) as a function of Pi and P 2 . In this 
way of plotting it is however rather difficult to see why the contribution with at 
least one scale above 1 GeV is as large as shown in Tab. HJ The quantity a^ L is 
much more suitable for this. The result for after integrating for this quantity is 
directly proportional to the volume under the surface as it is plotted in Figs. El 
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and [7] with a logarithmic scale for Pi and Pi. We have used the same scale for all 
three plots. What one finds is that the VMD one has much smaller contributions 
for Pi and Pi large but both MV and KN show a significant contribution even at 
fairly high values of (Pi, Pi). Also the contribution at these higher values of (Pi, Pi) 
is concentrated along the axis Pi = Pi. One also see by comparing Figs. [5] and [6] 
that the enhancement of the MV result over the KN result comes not from a very 
different shape but more a general increase over the entire region. The parts below 
0.1 GeV were not plotted, these are very similar for all three cases. A plot for the 
WZW case simply shows a constantly growing ridge along Pi = P% which produces 
then the log 2 A divergence. 

The figures before give an indication of which ranges of (Pi, Pi) are important. 
But what about the values of Q that are relevant. This will of course depend on the 
values of Pi and Pi . We show in Fig. [8] the value for a^ L< ^ along the line Pi = P%. 
Again, the contribution to is proportional to the volume under the surface as 
shown. This is shown for the MV and KN form factors in Figs. [8] and [9] respectively. 
One surprise for us was that while one can see that the tail towards larger values of 
Q is somewhat larger for the MV form factor than the KN one, it is much less than 
expected and only marginally visible in the plot. The main conclusion from this 
section is that the numerical difference between MV and KN comes from relatively 
low values of Q and moderate values of Pi and Pi. We have provided plots and 
numerics so that readers can draw their own conclusions. 
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5. Comparison 

Let us now try to compare the different results of the three calculations in Refs. 
[i [H EE CDS and OH In Tab. the results to leading order in 1/N C are shown. 
The quark loop is of the same order and has to be added to get the full hadronic 



I 



I 



2, 2008 8:43 WSPC/INSTRUCTION FILE gm2review 



The Hadronic Light-by- Light Contribution to the Muon g — 2: Where Do We Stand? 13 





light-by-light while the model used in Ref. Q~7]is saturated just by exchanges. In the 
GBE the effect of the new OPE in Ref. 17 is a little larger than the quark loop 
contributions of Refs. 181 and 1 131 but compatible within one sigma. This contribution 
has been discussed in more detail in the previous section. The new OPE in Ref. [T7| 
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Table 7. Full hadronic light-by-light contribution to 
a M at 0(N C ). The difference between the two re- 
sults of Refs. [8] and 1131 is the contribution of the 
scalar exchange —(0.7 ± 0.1) X 10~ 10 . This contri- 
bution is not included in Refs. 1101 1141 and 1171 



Hadronic light-by-light at 0(N C ) 


10 10 X a M 


Nonet Symmetry + Scalar Exchange ^ ^ 
Nonet Symmetry 
Nonet Symmetry l 10 l 14 l 
New OPE and Nonet Symmetry l^ZI 
New OPE and Ideal Mixing El 


10.2 ± 1.9 
10.9 ± 1.9 
9.4 ± 1.6 
12.1 ± 1.0 
13.6 ± 1.5 



similarly increases the axial-vector exchange with nonet symmetry from 0.3 xl0~ 10 
to 0.7 x 10~ 10 One thus sees a reasonable agreement in the comparison of the 0(N C ) 
results of Refs. 151 HHl ITDl ITU and H~71 when using the same mass mixing for the axial- 
vectors, namely, (10.9± 1.9, 9.4 ± 1.6, 12.1 ±1.0). 

The final differences are due to the additional increase of 1.5xl0~ 10 from the 
ideal mixing in the axial vector exchange in Ref. [17] and the scalar exchange of 
-0.7xl0- 10 in Refs. 8 and|H 

Let us now see what the different predictions at NLO in 1/N C are. In Ref. H~7l 
the authors studied the chiral expansion of the charged pion loop using the HGS 
model used in Refs. [TO] and [13] This model is known not to give the correct QCD 
high energy behavior in some two-point functions, in particular it does not fulfill 
Weinberg Sum Rules, see e.g. Ref.[8j Within this model, Ref. 17 showed that there is 
a large cancellation between the first three terms of an expansion of the charged pion 
loop contribution in powers of {m^/Mp) 2 . It is not clear how one should interpret 
this. In Refs.[8jsome studies of the cut-off dependence of this contribution were done 
and the bulk of their final number came from fairly low energies which should be 
less model dependent. However, it is clear that there is a large model dependence 
in the NLO in 1/N C contributions. But simply taking it to be (Oil) x 10~ 10 as in 
Ref. [XT1 is rather drastic and certainly has an underestimated error. The argument 
of very large higher order corrections when expanded in CHPT orders which was 
used against this contribution in Ref. 17] also applies to the ir° exchange as can be 
seen from Tab. [6] by comparing the WZW column with the others. 

Let us now compare the results for the full hadronic light-by-light contribution 
to a M when summing all contributions. The final result quoted in Refs. 18] ITS" ] ITO ] fT4l 
andfTTlcan be found in Tab. [8] The apparent agreement between Ref s . 18] [T3l and 1T0] 
[2] final number is hiding non- negligible differences which numerically compensate 
to a large extent. There are differences in the quark loop and charged pion and 
Kaon loops and Refs. [1111141 do not include the scalar exchange. 

Comparing the results of Refs . [51 [TBI and 1T71 we have seen several differences of 
order 1.5 x 10 -10 , differences which are not related to the one induced by the new 
short-distance constraint introduced in Ref. [T71 These differences are numerically 



I 



I 
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Table 8. Results for the full hadronic light-by-light 
contribution to a u . 



Full Hadronic Light-by-Light 


10 10 X a M 


Bijnens, Pallante and Prades ^ ^ 
Hayakawa and Kinoshita 
Melnikov and Vainshtein 


8.3 ± 3.2 
8.9 ± 1.7 
13.6 ± 2.5 



of the same order or smaller than the uncertainty quoted in Refs. [8] EH but tend to 
add up making the total difference large as follows: The different axial-vector mass 
mixing account for —1.5 x 10~ 10 , the absence of scalar exchange in Ref . [171 accounts 
for —0.7 x 10~ 10 and the absence of the NLO in 1/N C charged pion and Kaon 
loops contribution in Ref. [T7| accounts for —1.9 x 10~ 10 . These model dependent 
differences add up to —4.1 x 10~ 10 out of the final —5.3 x 10 -10 difference between 
the results in Refs. [13] and [TT] In addition we have shown from which regions 
in momentum the main contribution originates. Clearly, the new OPE constraint 
found in Ref. [17] alone does not account for the large final numerical difference with 
respect to Refs. [13] as a reading of it seems to suggest. 

6. Conclusions 

At present, the only possible conclusion is that the situation of the hadronic light- 
by-light contribution to a M is unsatisfactory. However, looking into the various calcu- 
lations one finds a numerical agreement within roughly one sigma when comparing 
the <D(N C ) results found in Refs. EHHIIH [Hand El see Tab. [3 A new full 0{N C ) 
calculation studying the full correlator with the large N c techniques developed in 
Refs. l22l and l23l and references therein, seems feasible and definitely desirable. 

At NLO in 1/N C , one needs to control both Goldstone and non-Goldstone boson 
loop contributions. The high model dependence of the Goldstone boson loop is 
clearly visible in the different results of Refs . 18] [13] and ITO] [141 and discussed in Refs. 
[8j and [TTJ For non-Goldstone boson loops, little is known on how to consistently 
treat them, a recent attempt in another context is Ref. [2H 

In the meanwhile, we propose as an educated guess for the total IiLBlE) 

a M = (11 ±4) x 10~ 10 . (11) 

We believe that, that this number and error capture our present understanding of 
the hLBL contribution to et M . This number can be reached using several different 
arguments: the new short-distance constraint found in Ref. 1171 and the ideal mixing 
for the axial- vector exchange should lead to some increase of the results of Refs. 
[TBI and [TOl [HI the scalar exchange and the pion and Kaon loops are expected to 

e This educated guess agrees with the one presented by Eduardo de Rafael 2 and ourselves E] a t 
the "Final Euridice Meeting" in Kazimierz, August 2006 and by one of us (JB) at the "DESY 
Theory Workshop" in Hamburg, September 2005. 
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lead to some decrease of the result of Ref. [17] one can also average the leading in 
1/N C results (three middle results of Tab. [7]) which turn out to be within one sigma. 
The final error remains a guess but the error in is chosen to include all the 
known uncertainties. 
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